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We study the QCD phase diagram, in particular we study the critical 
points of the two main QCD phase transitions, confinement and chiral sym- 
metry breaking. Confinement drives chiral symmetry breaking, and, due 
to the finite quark mass, at small density both transitions are a crossover, 
while they are a first or second order phase transition in large density. 
We study the QCD phase diagram with a quark potential model including 
both confinement and chiral symmetry. This formalism, in the Coulomb 
gauge hamiltonian formalism of QCD, is presently the only one able to 
microscopically include both a quark-antiquark confining potential and a 
vacuum condensate of quark-antiquark pairs. This model is able to address 
all the excited hadrons, and chiral symmetry breaking, at the same token. 
Our order parameters are the Polyakov loop and the quark mass gap. The 
confining potential is extracted from the Lattice QCD data of the Bielefeld 
group. We address how the quark masses affect the critical point location 
in the phase diagram. 

PACS numbers: PACS numbers come here 



1. Introduction 

Our main motivation is to contribute to understand the QCD phase 
diagram [1], for finite T and to be studied at LHC, RHIC and FAIR. 
Moreover, our formalism, in the Coulomb gauge hamiltonian formalism of 
QCD, is presently the only one able to microscopically include both a quark- 
antiquark confining potential and a vacuum condensate of quark-antiquark 
pairs. This model is able to address excited hadrons and chiral symmetry at 
the same token, and we recently suggested that the infrared enhancement of 
the quark mass can be observed in the excited baryon spectrum at CBELSA 
and at JLAB [2j i3j. Thus the present work, not only addresses the QCD 



* Presented at Excited QCD 2010, Tatranska Lomnica; work partly done with Tim 
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Fig. 1. (right) A sketch of the QCD Phase Diagram, according to the cohaboration 
CBM at FAIR [I], it is usuahy assumed that the critical point for deconfinement 
coincides with the critical point for chiral symmetry restoration, (left) First calcu- 
lation of excited baryons with a chiral invariant quark model. 



phase diagram, but it also constitutes the first step to allow us in the future 
to compute the spectrum of any hadron, extending, say the Fig. [1] (left) 
computed in reference ^ to finite T . 

Here we address the finite temperature string tension, the quark mass 
gap for a finite current quark mass and temperature, and the deconfinement 
and chiral restoration crossovers. We conclude on the separation of the 
critical point for for chiral symmetry restoration from the critical point for 
deconfinement . 

2. Fits for the finite T string tension from tlie Lattice QCD 

energy 

At vanishing temperature T = 0, the confinement, modelled by a string, 
is dominant at moderate distances, 

V(.r) c^-!^+Vo + ar . (1) 
12r 

At short distances we have the Luscher or Nambu-Gotto Coulomb due to the 
string vibration + the OGE coulomb, however the Coulomb is not important 
for chiral symmetry breaking. At finite temperature the string tension (7(T) 
should also dominate chiral symmetry breaking, and thus one of our crucial 
steps here is the fit of the string tension o"(r) from the Lattice QCD data 
of the Bielefeld Lattice QCD group, 01 EJ El [8] . 

The Polyakov loop is a gluonic path, closed in the imaginary time 
(proportional to the inverse temperature T~^) direction in QCD discretized 
in a periodic boundary euclidian Lattice. It measures the free energy F of 
one or more static quarks, 

P(0) = Ne-^^/^ , P«(0)P^(r) = Ne-^'"i'^'''>/^ . (2) 
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Fig. 2. (top left) The Bielefeld free Fi energy at T < T^. (top right) Detail of 
the string tension fit in the case of T = Q.9ATc- (bottom left) The critical curve 
for ^ ^ as a function of y^, for T ~ Tc it behaves like a square root, (bottom 
right) comparing the magnetization critical curve with the string tension cr/fToj 
fitted from the long distance part of Fi, they are quite close. 

If we consider a single solitary quark in the universe, in the confining phase, 
his string will travel as far as needed to connect the quark to an antiquark, 
resulting in an infinite energy F. Thus the 1 quark Polyakov loop P is a 
frequently used order parameter for deconfinement. With the string tension 
cr(T) extracted from the qq pair of Polyakov loops we can also estimate 
the 1 quark Polyakov loop -P(O). At finite T, we use as thermodynamic 
potentials the free energy Fi and the internal energy C/i, computed in Lattice 
QCD with the Polyakov loop [H El El O [8] . They are related to the static 
potential V{r) = —fdr with Fl(r) = —fdr — SdT adequate for isothermic 
transformations. In Fig. [2] we extract the string tensions cr(T) from the free 
energy Fi{T) computed by the Bielefeld group, and we also include string 
tensions previously computed by the Bielefeld group [9]. 

We also find an ansatz for the string tension curve, among the order 
parameter curves of other physical systems related to confinement, i. e. 
in ferromagnetic materials, in the Ising model, in superconductors either 
in the BCS model or in the Ginzburg-Landau model, or in string models. 
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Fig. 3. (left) The mass gap m(0) solution of as a function of the quark current 
mass mo, in units of cr — 1. (right) Sketch of the effect of mo on the crossover 
versus phase transition of choral restoration at finite T. 



to suggest ansatze for the string tension curve. We find that the order 
parameter curve that best fits our string tension curve is the spontaneous 
magnetization of a ferromagnet [lOj, solution of the algebraic equation, 



M 



Tc M 



tanh . (3) 

In Fig. [2] we show the solution of Eq. [3] obtained with the fixed point 
expansion, and compare it with the string tensions computed from lattice 
QCD data. 



3. The mass gap equation with finite T and finite current quark 

mass mo. 

Now, the critical point occurs when the phase transition changes to a 
crossover, and the crossover in QCD is produced by the finite current quark 
mass mO, since it affects the order parameters P or a, and the mass gap m(0) 
or the quark condensate (qq). The mass gap equation at the ladder /rainbow 
truncation of Coulomb Gauge QCD in equal time reads, 

a r°° dk Ia{p-, k, n) m{k)p — Ib{p, k, fi) m{p)k 

10 



m{p) 
lA{p,k,n) 
lB{p,k,fi) 



mo + ^ 



27r 
pk 



y/k"^ + m{k)^ 
pk 



(4) 



{P 



k)^ + fi^ 
pk 



{p - k)' 



+ 



(p + A:)2 + ^2 

pk 1 (p 

ip + k)^ + ^i^ + 2 °^ {p + k)-^ + fi^ 



k)^ + fl^ 



The mass gap equation (UD for the running mass ■m{p) is a non- linear integral 
equation with a nasty cancellation of Infrared divergences |1H IT^ |T3] . We 
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Fig. 4. Sketches of the saturation of confinement (left), and of the corresponding 
crossover in the order parameter P polyakov loop (righ). 

devise a new method with a rational ansatz, and with relaxation, to get a 
maximum precision in the IR where the equation is nearly almost unstable. 
The solution m{p) is shown in Fig. [3] for a vanishing momentum p = 0. 

At finite T, one only has to change the string tension to the finite T 
string tension (j{T) [H], and also to replace an integral in ci; by a sum 
in Matsubara Frequencies. Both are equivalent to a reduction in the string 
tension, o" — >■ a* and thus all we have to do is to solve the mass gap equation 
in units of a* . The results are depicted in Fig. [3l Thus at vanishing mo 
we have a chiral symmetry phase transition, and at finite rriQ we have a 
crossover, that gets weaker and weaker when mo increases. This is also 
sketched in Fig. [3l 

4. Chiral symmetry and confinement crossovers w^ith a finite 

current quark mass 

In what concerns confinement, the linear confining quark- antiquark po- 
tential saturates when the string breaks at the threshold for the creation of 
a quark-antiquark pair. Thus the free energy -F(O) of a single static quark 
is not infinite, but is the energy of the string saturation, of the order of the 
mass of a meson i. e. of 2mo. For the Polyakov loop we get, 

P(0) ~ Are-2™o/^ . (5) 

Thus at infinite mo we have a confining phase transition, while at finite mg 
we have a crossover, that gets weaker and weaker when mg decreases. This 
is sketched in Fig. HI 

Since the finite current quark mass affects in opposite ways the crossover 
for confinement and the one for chiral symmetry, we conjecture that at 
finite T and n there are not only one but two critical points (a point where 
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a crossover separates from a phase transition). Since for the light u and 
d quarks the current mass mo is small, we expect the crossover for chiral 
symmetry restoration critical to be closer to the /x = vertical axis, and 
the crossover for deconfinement to go deeper into the finite /i region of the 
critical curve in the QCD phase diagram depicted in Fig. [TJ 

5. Outlook 

We soon plan to complete the part of this work which is only sketched 
here, i. e. to compute the crossover curves for the chiral symmetry restora- 
tion and for the deconfinement at finite T. This will also require the study 
of light hadrons at finite temperature T. Thus we also plan to address the 
excited hadron spectrum at finite T, continuing the work initiated with Tim 
Van Cauteren, Marco Cardoso, Nuno Cardoso and Felipe Llanes-Estrada. 
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